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According to recent experimental data at GSI, the rate of the number of daughter ions ^''"Ce^*^, produced by 
the nuclear K-shell electron capture (EC) decay of the H-like ion i'*0pj.58+^ jg modulated in time with a period 
Tec = 7.06(8) sec and an amplitude uec ~ 0.18(3). We show that this phenomenon can be explained by neutrino 
mass differences and derive a value for the difference of squared masses Am2i = mi — nii — 2.22(3) X 10-* eV^ 

PACS: 12.15.Ff, 13.15.+g, 23.40.Bw, 26.65.+t 



Introduction 

The experimental investigation of the K-shell 
electron capture (EC) and f3~^ decays of the H- 
like i40pj.58+ ^j^g He-like i40pj.57+ ^^^^^ 

been recently carried out in the Experimental 
Storage Ring (ESR) at GSI in Darmstadt [1]. The 
obtained results showed a dependence of the weak 
decay rates on the electron structure of the heavy 
ions. As has been shown in |2j, the experimen- 
tal data on the ratios of the weak EC and f3~^- 
decay rates of the H-like i40pj.58+ 
like I40pr57+ jQ^g^ obtained in GSI [T], can be 
described within standard theory of weak interac- 
tions of heavy ions and massless Dirac neutrinos 
[3] with an accuracy better than 3%. 

However, a very recent measurement of the 
time-dependence of the rate of the number of 
daughter ions ^40Ce^^"'" from the i?C-decay of the 
H-like i40pr58+ ion, i.e. I40pr58+ ^ ^^^^^58+ _^ 

Ve, showed a time-modulation of the exponential 
decay with a period Tec = 7.06(8) s and an am- 
plitude uec = 0.18(3) 4JI. Since the rate of the 
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perimental value of the period of the time— modulation 
is Tec = 7.10(22) s with the amplitude of the time- 
modulation asc = 0.23(4) 13]. 



number of daughter ions is defined by 



it) 



dt 



(1) 



where A_ec(0 is the iJC-decay rate and Nra{t) is 
the number of mother ions i40pj.58-i-^ ^j^^ time- 
modulation of dNf'-^{t)/dt implies a periodic 
time-dependence of the i?C-decay rate XEc{t) 

m 



XEc{t) = Xec \ 1 + clec cos I 



27ri 
Tec 



^Ec)}, (2) 



where aEc, Tec and 4>ec ai'e the amplitude, pe- 
riod and phase of the time-dependent term [3]. 

Nowadays the existence of massive neutrinos, 
neutrino-flavour mixing and neutrino oscillations 
is well established experimentally and elaborated 
theoretically [5]. However, these phenomena con- 
cerned mainly the propagation of solar neutrinos 
and reactor anti-neutrinos in space and time com- 
bined with neutrino oscillations [5 . Wc show that 
the observed time-modulation of the rate of the 
number of daughter ions dNjj^^ {t) / dt in the EC- 
decay of the H-hke i40pj.58+ -^^^^ caused 
by neutrino oscillations but can be explained by 
a kind of quantum beats [B] due to mass differ- 
ences of neutrino mass-eigenstates [5]. This can 
provide a new method for studying of massive 
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neutrino mixing, neutrino mass differences and 
neutrino vacuum polarisation. 

Amplitudes of E'C decays of H like heavy 
ions 

The Hamilton operator Hw(t) of the weak in- 
teractions, responsible for the EC and /3+ decays 
of the H-like heavy ions, can be taken in the stan- 
dard form [5] but accounting for the neutrino- 
flavour mixing [5] (see also l7j-[9]). This gives 
^w{t) — J^j^eji^w (t), where Uej are matrix 
elements of the neutrino mixing matrix U [S]. 
The Hamilton operator of the weak interactions 
B.\^{t) is defined by 



x[V5.,(x)7M(l-7')V'e-(a;)], 



(3) 

with standard notations [2|3j. 

The massive neutrinos Vj in the final state of 
the -EC-decay m — > d + where m and d are 
the mother and daughter ions, are indistinguish- 
able in principle, since the electron is entangled 
with the electron neutrino only ^ , which is the 
superposition of the neutrino mass-eigenstates 
— 15]. Such an indistinguisha- 

bility of massive neutrinos Vj requires to take the 
amplitude A{m d+i'e){t) of the E'C-decay as a 
coherent sum of the amplitudes A{m d + i'j){t) 
of the m d -\- Vj decays, describing three alter- 
native ways of the evolution of the initial system, 
i.e. the mother ion m into the daughter ion d. 
According to Feynman jlQj, the probability am- 
plitude of the evolution of the quantum system 
should be taken in the form of a coherent sum 
of the probability amplitudes of every alternative 
way of the evolution (sec also llj) in contrast to 
the recent assertion by Giunti |12j and Kienert et 
al. [13] . As a result we get 

A{fn -> d + Ue){t) = ^ UejA{m d + i^j){t), (4) 



where the coefficients Uej take into account that 
the electron couples to the electron neutrino only. 
In the standard time-dependent perturbation 



theory [Tl] the amplitudes A{t 



given by 

A{m ~> d + Vj){t) 



-I I dre^'' {diqy,ik,)\Hll\T)\miO)), (5) 

' — oo 

where q and kj are 3-momenta of the daughter 
ion and the neutrino mass-eigenstate Vj, respec- 
tively. For the regularization of the integral over 
time we use the e — > regularization procedure. 

For the description of the time modulated inter- 
ference term in the E'C-decay rate of the H-like 
heavy ion Eq.([2]) we assume a non-conservation 
of 3-momenta of neutrino mass-eigenstates in or- 
der to deal with differences of 3-momenta of mas- 
sive neutrinos playing important role in our anal- 
ysis of the interference term. Non-conservation of 
3-momenta of massive neutrinos can be caused 
by the uncertainties of momenta of the bound 
electron, proton and neutron in the elementary 
e~ -\- p n + Vj transition of the i?C-decay of 
the H-like heavy ion (see also [15j). 

A quantum mechanical description of the EC- 
decays of the H-like heavy ions with different 
3-momenta of neutrino mass-eigenstates can be 
carried out using the wave functions of massive 
neutrinos in the form of the Gaussian wave pack- 
ets 1161 



V'.,(r,t)-(27r52)3/2 



d^k 



(27r)3 

xe^^-^^-*^^(^)^^.^.(^,a.^.), 

where (5 is a spatial spread of the massive neutrino 
Vj, kj is the neutrino momentum and Ej{k) = 



(6) 



fc 2 + m?- is the energy of a plane wave with the 

momentum fc, u^. (fc, a^.) is the Dirac bispinor of 
the massive neutrino Vj. In the limit i5 ^ oo the 
wave function ([6]) reduces to a plane wave. 

Following [2], we obtain the amplitude of the 
i?C-decay as a function of time t 



A{m -> d -f Ve){t) = -V3y^2M^2E4q) 



iA2/;r , l:.^2 pii^Ej{q)~ie)t 



defined in the rest frame of the mother ion, are 



X e 2 



^^,SHq-+k,ye_^ 



AE,{q) 
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where AEj{q) = Ed{q) + Ej{q) - M„ is the 
energy difference of the final and initial state, 
q is the S-momentum of the daughter ion d, 
Ed{q) = \/q"^ + and are the energies 
of the daughter and mother ions, respectively, 
and Ej{q) is the energy of the neutrino mass- 
eigenstate vj with momentum q, A^gt is the nu- 
clear matrix element of the Gamow-Teller tran- 
sition m ^ d and ("01^ ) is the wave function of 
the bound electron in the H-like heavy ion m, 
averaged over the nuclear density [|2] . 

The rate of the neutrino spectrum of the EC- 
decay as a function of time is defined by 



dt 

1 



1 



d\ 



2F+ldt 



2M„y i2n)^2Ed{q) 



(8) 



The integration over q can be carried out if the 
width S of the wave packets of the wave functions 
of neutrino mass-eigenstates is sufficiently large, 
so that the Gaussian function e ~ * (q+p) jg lo- 
calised in the vicinity of (f ~ ~p, where p — kj 

for the diagonal term and p = (ki + kj)/2 = k^^-' 
for the interference term. The result of the inte- 
gration over q is given by two terms 



dN^^ jt) _ dN!)l\t) dNl">{t) 
dt dt dt ' 

where we have denoted the diagonal term 



(2), 



(9) 



dNi'-Xt) 



dt 



2F 



2e 



{AE,{k,W+e^ 



e^^ (10) 



and the interference term 



EMf)E,{m^) 



2e 



2e 



(+)^ 



12 I c2 



,et 



X cos 



(eS'i^^)-e,{}^^^) 



(11) 



Here k\j^ = {k, - kj)/2 and k\f^ = {k, + kj)/2 
are the difference and averaged neutrino momenta 

and Ej{k\f) = ^{k'^f)^+m]. The former is 
an approximate relation between the energy of 
the massive neutrino Vj and a momentum kl^\ 
since 3~momenta of massive neutrinos are not 
conserved and massive neutrinos are off-shell. As 
a result such a relation cannot be used for the 
analysis of energy differences Ei {k^^^ ) — Ej (k^^^ ) , 
which are sensitive to the off-shell and on-shell 
states of massive neutrinos. We show this below. 

The EC— decay rate Xec of the H like ions 

The diagonal part of the rate of the neutrino 
spectrum dNll\t)/dt defines the meanvalue of 
the £'C-decay rate = Xec- Taking the limit 
£ in Eq.Hni) we get 

r^'\t) : 



dm 



dt 



2F + 1 



\Mot\' {7.6^/^ 



X 2_^\Ue,\'E,{k,)2ndiAE,ik,)), 



(12) 



where the (5-functions d{AEj(kj)) with the argu- 
ments AEj{kj) = Ed{kj) + Ej{kj) — Mm describe 
the conservation of energy in the i?C-decay chan- 
nels m — > d + Vj. The solution of the equations 
AEj{kj) = gives the energies of neutrino mass- 
eigenstates shown in Fig. 1 with the energy dif- 



ferences Ei{ki) — Ej{kj) 



Am^j/2M„ 



For the calculation of the i?C-decay rate Xec 
we can set neutrino masses zero and use the uni- 
tarity of the [/-matrix ^ ■ |J7ejP = 1- This gives 



d^k 



dN!.^\t) 



{2TTf2E^ (7rJ2)3/2 



dt 



(13) 



where E^ — \k \ and (TrJ^)^/^ is related to the nor- 
malisation of the neutrino wave functiorQ. Hav- 

"^For sufficiently large S one gets f d?x i/it (r, (r, t) = 

(7r52)3/2e-*^(^'V^)'' 2_Bj. In the limit <5 ^ oo the 

Gaussian function behaves as (wS^)^/'^ e~ ^ ^ 2 — > 
(27r)'^ (5(^)((fj — kj) . This leads to the standard relativis- 
tic covariant normalisation of the neutrino wave function 
f cfix'^i(r,t)ip^ (r,t) = (2TTf2Ej5^^\qj - kj). For a 

finite S and qj = kj we get / d^xi/il {r,t)'<p^ (T,t) = 

(7r52)3/2 2£;^.. 
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Figure 1. Schematic diagram of the neutrino en- 
ergies Ej{kj) = {M^ - Mj + m2)/2M„ ~ Qh, 
where Qh = Afm — Md is the Q-value of the EC- 
decay. The neutrino energy differences are equal 
to Ei{ki) - E,{kj) = = Amy2M™. 



ing integrated over the neutrino phase volume we 
obtain [2] 

^EC^^^l\MoT\\4f)?^, (14) 

where Qh = 3.348(6) MeV for the ^JC-decay 
i40pr58+ ^ i40cg58+ ^ ^ The i;C-decay 

rate \ec has been calculated in [2] (see also [11]) 
within standard theory of weak interactions of 
heavy ions [3]. Since experimental value of the 
matrix element Uiz = sin ^13 e * "^^^ , where 5c p 
is a CP-violating phase [5], is very close to zero, 
we set = and below deal with two neu- 
trino mass-eigenstates only. The elements of the 
mixing matrix are set equal to Uei — cos 9i2 and 
= sin0i2 0. 

i?C— decays of H like heavy ions as analog 
of quantum beats of atomic transitions 

The time-modulation of the i?C-decays of the 
H-like heavy ions bears similarity with quantum 
beats of atomic transitions [B^, since in the EC- 
decays one deals with the transitions from the 
initial state \m) to the final state |fii^e), where 
the electron neutrino is the coherent state of two 
neutrino mass-eigenstates with the energy differ- 
ence equal to bJ2i = Arn^i/^Mm- The time dif- 
ferential detection of the daughter ions from the 
SC-decays in the GSI experiments with a time 
resolution ~ 1 s introduces an energy uncer- 



tainty 6 Ed 2TTh/Td = 4.14 X 10-15 gv. Thus, if 
6Ed > UJ21, following the analogy with quantum 
beats of atomic transitions [B] one should expect a 
periodic time-dependence of the i?C-decay rate 
with a frequency UJ21 and a period Tec ~ Mm- 

Analysis of the interference term of EC- 
decay rates of H like heavy ions 

The dependence of the interference term on the 
average momentum A:2i~^ = (^2 + ki)/2 can pro- 
duce an impression that the frequency of the pe- 
riodic time-dependence of the interference term 
is of order of $^21 — ^rn^il'^Qti o^^ly [12] 

cos{{E2{P+^)-E,{k!i'>))t)^ 
/ m\ — m\ \ 

'i52(4tVi?i(fc^t^)*' ^ 
=^cos(4^i) =cos(l]2i0, (15) 

where Ei{k^^^) ~ £^2(^21^^) - Qn and on-shell 
relations between the energies and momentum 
have been used [T7]. However, this result, being 
only partly correct, leads to the missing of the 
frequency 0^21 due to the use of on-shell relations 
between energies and momenta for the energy dif- 
ference. The existence of the frequency a'21 can 
be shown using non-conservation of neutrino 3- 
momenta, their differences and noticing that the 
behaviour of the interference term is governed by 
the Gaussian function e ^ 21 ) . 

The region of variation of 3-momenta of mas- 
sive neutrinos can be divided into two parts with 
the absolute values | /c2 1 > \ki and | fc2 1 < \ki\, re- 
spectively, and both cases lead to totally different 
modulation frequencies of the interference term. 

For 1^2! > 1^1, the energies i?i(fci), i?i(fc2i), 
£^2(^21) and £^2(^2) satisfy the inequality: 
Eiih) < Ei{k2i) < E2{k2i) < E2{k2) as it 
shown in Fig. 2. The deviations as indicated in 
Fig. 2 are of order of 0{\/5). Indeed, for suffi- 
ciently large (5, that has been already assumed 
above for the calculation the rate of neutrino 
spectrum Eq. pi}]) and Eq. pTj) . due to the func- 
tion e - ■'^^^^ai the region of the 3-momenta of 
massive neutrinos is constrained by /c2 ~ ^i- 
The allowed region of deviations is of order of 



Figure 2. For \k2\ > |fci| (left) the interfer- 
ence term has a periodic time-dependence with 
a frequency UJ21 due to the relation Ei{ki) < 
Ei{ki+^) < ^2(41^) < ^2(^2) (right) with 
kiV - (fc2 + fci)/2. 



Figure 3. For \k2\ < \ki\ (left) the interference 
term has a periodic time-dependence with a fre- 
quency fl2i > UJ21 due to the relation £'1(^21'^) < 
Ei{ki) < E2{k2) < E2{ki'^^) (right) with = 
(fc2 + fci)/2. 



1^1 ''1 ~ 1/^- ^'^^ these momenta the energies of 
neutrino mass-eigenstates obey the obvious rela- 
tion E2{k2) ~ Ei{ki) > |fc^7^| - 1/(5. Making 
the expansions of the energies of massive neutri- 
nos in powers of fcjj^ ^ and keeping only the first or- 
der contributions we get £^2(^:2^'') = £2(^2) — A£2 
and £i(A?^|^) = Ei{ki) + A£i, where we have 
denoted A£2 — fc2 • 47V'^2(^2) and A£i = 
ki-k2i'' / El {ki ) . Without loss of generality we can 
set A£2 ^ A£i ~ A£, where \ A£\ = 0(1/(5) and 
£2(^2) - Ei{ki) > \A£\ = 0{l/5). This shows 
that the deviations of £2(^21^'') and Ei{k^^) from 
£2(^2) and Ei{ki), respectively, arc of order of 
0(l/(5) and the interference term should have a 
periodic time dependence with a period UJ21- 

In the region \k2\ > |fci| of 3-momenta of mas- 
sive neutrinos the calculation of the periodic term 
runs as follows 

C08{{E2{kit^)-EMV))t)^ 

^ cos((£2(fc2) - Ei{ki))t) = cos(c^2ii), (16) 

where we have used the on-shell relations between 
neutrino energies and momenta only for the ex- 
pansions E2{kii^) = £2(^2) - A£ ^ E2{k2) - 
0{l/5) and £;i(4t^) Ei{ki) + A£ = Ei{ki) + 
0{l/5) but not for the energy difference. This 



gives E2{kfi^) - Ei{ki+^) = £2(^2) - Eiih) - 
2A£ = uj2\ — 0{2/5). As a result the validity 
of the interference term to have a frequency u>2i 
should be confirmed by the inequality ^21 2> 2/(5, 
which can be always satisfied for the sufficiently 
large 5. 

In turn, for \k2\ < |fci| the energies i?i(fc2i^''), 
Ei{ki), £2(^2) and E2{k2i^) satisfy the inequal- 
ity: Ei{ki+^) < £i(fci) < £2(^2) < E2{kit'^) as 
it shown in Fig. 3. In this case a periodic time- 
dependence of the interference term can be de- 
fined by the frequency, which is greater than 1021- 

In order to take into account both frequencies 
of the periodic time-dependence of the interfer- 
ence term we introduce the function 

p(fc2,fcl) = 0(|fc2|-|fcl|)p'(fc2,fcl) 

+ 0(|fcl|-|fc2|)p"(fc2,fcl), (17) 

having a meaning of the probability density for 
neutrino mass-eigenstates z/2 and Ui to have 3- 
momenta ^2 and ki, respectively. The Heavi- 
side step functions 9{\k2\ — |fci|) and 9{\ki\ — \k2\) 
take into account two possibilities 1^2! > l^il and 
1^2! < 1^1 1 for the 3-momenta of neutrino mass- 
eigenstates shown in Fig. 2 and Fig. 3. 

The contribution of the interference term 
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'^BC (^) -EC-decay rate we define as follows 



X p(fc2,fci) 



(7r52)3/2 



dt 

d^k2 



(18) 



(2^)32i;i(fci) (27r)32i;2(fc2) 
Substituting Eq.^ into Eq.([T8l) we get 

^^Ekt) = ^Ecit) + ^Ecitl (19) 

where {t) and (t) are determined by the 
contributions of the regions |fc2| > l^il and 1^21 < 
I fci I , respectively. They are given by the following 
momentum integrals 



Ag^(t) - mil 



2F 



\Mgt\'\{4T)\' 



i sin 26*12 



(2^)32£;i(fci) (2^)32i;2(fc2: 



■^^(|fc2|-|fcl|) 



X p'(fc2, fci) e 

2e 2e 



(AE^ikii^))^ + £2 (A£;i(fc^+)))2 + £2 



(20) 



and 

Ag^'(i)-lim^ 



.^c: y^J - I-^gtI^ 2 ^^"2^12 



d3fci 



(27r)32£;i(fci) {2TT)^2E2ik2 



■^(|fci|-|fc2|) 



xp"(fc2,fci)e-^'(^27^)^ V^2(4l^)i?i(4l^) 
2e 2e 



X cos 



E2(ki+^^ 



^21 



i?l(fc(|)^ 



(21) 



According to Figs. 2 and 3, the frequencies in two 
integrals Eqs. ([20|) and (|2T|) should be equal to 1^21 
and ri2i, respectively. 

Making expansions in powers 

of 1^7^ I - i/sE 

and taking the limit e we reduce the r.h.s. 



^Formally this leads to the replacements Ai?2(fc2i'') 
AE2{k2), AEiik^^V) 



and -Ei(4t^) 



Ei{ki). 



of Eqs.dlOD and ^ to the form 



Ag^(t) 



2F+ 1 
d^ki 



\MoT\'\{4T)\'l,sm20,2 



'GTl 
d3fc2 



■9{\k2\-\k,\) 



(2^)32£;i(fci) {27T)^2E2{k2 
X p\k2, fci) e -^^i^ii? J E2{k2)Ei{h) 



2tt S{AE2{k2)) + 2nd{AEi{ki)) 

E2{k2) - Ei{ki)]t 



(22) 



and 



{t) 



2F+ 1 
d^ki 



sin 2f?i2 



d^k2 



(2^)32£;i(fci) {2^Y2E2{k2) 
;2^r(-)\2 



^(|fcl|-|fc2|) 



p"(^2,^i)e-^'(^2i )y E2ik2)Ei{h) 



2tt 5{AE2{k2)) + 27r(5(A£;i(A:i)) 



X cos 



E2{k2) + Ei{ki 



(23) 



By virtue of the J-functions (5(Ai?i(fci)) and 
5{AE2{k2)) with arguments A£;i(A:i) = Ed[ki) + 
Ei{ki) - Mm and Ai;2(fc2) = Edfy + Ei{k2) - 
Mm, respectively, providing energy conservation 
in the i?C-decay channels m — + d + vi and 
m ^ d + V2, and the function e^'^^'-'^^i the 
frequencies of the periodic terms can be replaced 
by CJ21 and O21, respectively. As a result the con- 
tribution of the interference term reads 



aEc cos{uj2it) + ttEc cos(r22ii), (24) 



Abc 



where Xec is defined by Eq. lfM]) . The amplitudes 
a EC and uec are given by the integrals 



O'EC — sin 201 



27r2 



d3A;i 



d^k2 



in J (27r)32i;i (fci) (271)32^2(^2) 
p'(fc2, fci) e -'^'(^21 y£^(^)^^ 

,5(A£;2(fc2)) + 6{AEi{ki))\ e{\k2\ - \ki\) (25) 
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and 

. 27r2 f d^h d'k2 

UEC = S\n20i2 -T-^ / — 

Q\ J {27T)^2Ei{ki) {2nf2E2{k2) 

X p"(k2, fci) e -^^&^^^E2{k2)Ei{h) 

X 5(AS2(fc2)) + <5(A£;i(^i))] |fc2|). (26) 

The amplitudes of the periodic functions, depend- 
ing on the mixing angle 9i2 and the parameter 
S, are considered as empirical parameters deter- 
mined by the experiment. 

We would like to emphasize that due to non- 
conservation of 3-momenta the massive neutrinos 
are off-shell. This means that on-shell relations 
Ej{p) — {p^ + ■m^Y^'^ between energies Ej{p) 

and 3-momenta p — kj or p — fc|^^ are only 
approximate. This results in some constraints 
on the use of on-shell relations between energies 
and momenta. Indeed, for the correct calcula- 
tion of the periodic time-dependence of the EC- 
decay rate in the momentum integrals Eqs. ipD)) - 
(PS)) one cannot use on-shell relations between 
energies and momenta of massive neutrinos for 
the analysis of the differences of neutrino ener- 
gies and momenta due to substantial sensitivity 
of these differences to the off-shell and on-shell 
states. As we have shown above (see Eq. p3)) and 
a discussion below) , the use of the on-shell rela- 
tions between the energies E2ik^i^^) and ^i(4l^) 
and the momentum fcji"^ leads to the missing of 
the frequency 0^21 (see [TH). 

Time dependent i?C decay rates of H like 
heavy ions 

Taking into account the contribution of the in- 
terference term, the total i?C-decay rate we get 
in the form 



A 



EC 



^ 1 + UEC cos{uj2it) + CLEC cos(r22it). (27) 



The i?C-decay rate Eq. (P7)) contains two periodic 
terms. In the laboratory frame the periods are 
equal to 

27r7 _ 27r7 

W21 "21 

where 7 = 1.43 is the Lorentz factor of the 
motion of mother ions in the laboratory frame 



Tec — 



(28) 



[5]. We identify Tec = 2717/^21 ~ A^m with 
the experimental period of the time-modulation 
Tec — 7.06(8) s, because it is shown experimen- 
tally that the period of the time-modulation does 
not depend on the Q-value of the weak transi- 
tion. For the £;C-decay of the H-like ^^^VTCi^°+ 
ion with Qh — 4827 keV the period of modula- 
tion Tec = 7.10(22) s is equal to Tec = 7.06(8) s 
of the H-like i40pj,58-i- ^^^^ within the experimen- 
tal error bars. Since the mass difference of the 
H-hke ions "2pni60+ and i40pj.58+ g^^jj 

expect only small differences in the periods of the 
time-modulation as shown by the experiment. 

Thus, for Tec = 7.06(8) s we get (Am|i)Gsi = 
2.22(3) X lO^-^eV^ (see also The value 

(Ato^i)gsi = 2.22(3) X lO'^'eV^ is by a factor of 
2.75 larger than (Am2i)KamLAND = O.SOig];^ x 
10^* eV^, obtained as a best-fit of the global 
analysis of the solar-neutrino and KamLAND ex- 
perimental data [5|9| . 

For (Aml^Gsi = 2.22(3) x lO^^eV^ we get 



1.8 X 10 s. These fast oscillations, aver- 



aged over the experimental time resolution AT = 
0.32 s !4i, are not observable. As a result the 
time-dependent i?C-decay rate Eq. (l771) is given 
only by the low frequency term 

A_Bc(0 = ^EC (1 + CLEC cos{LUEct)), (29) 

where luec = ^21 — '^.'n: /Tec and the amplitude 
aEc is defined by Eq. ([25l) . 



Discussion and summary 

We have shown that the experimental data on 
the time-modulation of the rate of the number 
of daughter ions ^^°Ce^'*^, observed in the EC- 
decay of the H-like ion i40pj.58+ -^^ 
plained by the neutrino mass differences. How- 
ever, the difference of squared neutrino masses 
(Amii)Gsi = 2.22(3) x lO^-^eV^, derived from 
the period Tec = 7.06(8) s of the rate of the num- 
ber of daughter ions ^^°Ce^^^, is 2.75 times larger 
than the value measured by KamLAND [5] . A so- 
lution of this problem in terms of neutrino mass 
corrections, induced by the interaction of mas- 
sive neutrinos with a strong Coulomb field of the 
daughter ion through virtual i~W~^ pair creation, 
is proposed in [15^. 

In summary we argue that the mechanism of 
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the time-dependence of the -EC-decay rates of 
the H-Hke heavy ions bears similarity with quan- 
tum beats of atomic transitions [5], as the EC- 
decays are the transitions of the initial state |m) 
into the final state {di^e), where the electron neu- 
trino is the coherent state of two massive neutri- 
nos with energy difference uj2i- A sophisticated 
calculation of the interference term shows also the 
existence of a periodic time-dependence with a 
frequency Q21 3> uj2i- In order to take into ac- 
count the contributions of two possible frequen- 
cies of the periodic time-dependence of the inter- 
ference term we have introduced the probability 
density p(fc2, fci) for neutrino mass-eigenstates 1^2 
and vi to get 3-momenta k2 and ki , respectively. 
The calculation of the function p{k2, ki) is rather 
hard problem. We can only argue that in the 
massless limit it should be equal to 

p(fc2 , fci ) = 2 ^E,iki)E,{k2){2TrfS^^^ (fe - fci ) , 

which is required by the necessity to retain the 
meanvalue of the EC-decay rate Xec- 

The necessary condition for the appearance of 
the interference term in the EC-decay rate is the 
overlap of the energy levels of neutrino mass- 
eigenstates. Non-conservation of 3-momenta of 
neutrino mass-eigenstates provides a possibility 
to place the interference term of the EC-decay 
rate in the energy region of the diagonal term, 
defining the meanvalue of the EC-decay rate 
Xec- a validity of such a transformation is sup- 
ported by the constraint W21 = Am2i/2Mm S> 
2/(5, which makes impossible the limit Am^i — > 
in the argument of the interference term, where 
the terms of order 0{l/5) are dropped with re- 
spect to Am2i/4Mm. This means that in the 
interference term of the EC-decay rate Eq. (|29|) 
one cannot set Ani^i — s- in order to reduce the 
time-dependent EC-decay rate to the EC-decay 
rate Xec given by Eg. p^ . A correct reduction of 
XEcit) to Xec can be carried out only by means 
of the averaging over time (XEcit)) = Xec- In 
spite of momentum non-conservation, energy of 
the EC-decay is conserved in our approach. This 
is required by Fermi Golden Rule 14J. 
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